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Abstract 



In supersymmetric theories of nature the Higgsino fermionic superpartner of the Higgs 
boson can arise as the lightest standard model superpartner depending on the couplings 
between the Higgs and supersymmetry breaking sectors. In this letter the production 
and decay of Higgsino pairs to the Goldstone fermion of supersymmetry breaking and the 
Higgs boson, h, or gauge bosons, Z or 7 are considered. Relatively clean di-boson final 
states, hh, hj, hZ, Zj, or ZZ, with a large amount of missing energy result. The latter 
channels provide novel discovery modes for supersymmetry at high energy colliders since 
events with Z bosons are generally rejected in supersymmetry searches. In addition, 
final states with real Higgs bosons can potentially provide efficient channels to discover 
and study a Higgs signal at the Fermilab Tevatron Run II. 



1 Introduction 



Supersymmetry (SUSY) provides perhaps the best motivated extension of the Standard Model. 
Spontaneous SUSY breaking leads naturally to radiative electroweak symmetry breaking with 
masses of order the electroweak scale for the superpartners of the Standard Model (SM) parti- 
cles. If the messenger interactions which couple the SM superpartners to the SUSY breaking 
sector are stronger than gravity, the lightest supersymmetric particle (LSP) is the Goldstone 
fermion of supersymmetry breaking, the Goldstino G. The next to lightest supersymmetric 
particle (NLSP) is generally the lightest SM superpartner. If the intrinsic scale of supersym- 
metry breaking is below ~ 10 3 TeV the NLSP can decay to its SM partner and the Goldstino 
on laboratory length scales [|TJ. This has an important impact on experimental SUSY sig- 
natures at high energy colliders. Since superpartners are generally produced in pairs, these 
decays give rise to final states with two hard partons and missing energy {$t) carried by the 
Goldstino pair, and with possibly other partons in the final state from cascade decays to the 
NLSP §0 0. 

The identity of the NLSP determines the type of final states which arise from decay to the 
Goldstino ||. A neutralino NLSP, Xi> which is gaugino-like, can decay by Xi ~~ > jG, leading 
to final states with 77$r- A slepton NLSP, £, can decay by £ £G, giving ££$t final states. 
In this letter we consider in detail the possibility of a fermionic Higgsino-like neutralino NLSP. 
Because it is the superpartner of the Higgs boson, h, a Higgsino NLSP can decay by Xi hG. 
In addition, since the longitudinal component of the Z boson mixes with the Goldstone mode 
of the Higgs field, Xi ~^ ZG can also result. Because of a strong phase space suppression of the 
h and Z final states near threshold, decay to a photon can also be important for Higgsinos not 
too much heavier than the Z boson. Pair production of Higgsinos which decay to Goldstinos 
can then give rise to the di-boson final states (hh, hj, hZ, Z7, ZZ)]pr @]- 

Di-boson signatures which include Higgs and Z bosons and I^t are quite novel discovery 
modes for supersymmetry in the mass range accessible to the current generation of high energy 
collider experiments. In conventional SUSY signatures, in which the lightest neutralino, x\i is 
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assumed to escape the detector without decay to the Goldstino, the mass splittings between 
supersymmetric particles required in order for h or Z to arise in a cascade decay, typically 
imply the superpartners are too heavy to be produced in sufficient numbers at present colliders. 
For this reason events with reconstructed Z bosons are in fact generally rejected in present 
SUSY searches. However, since the Goldstino is essentially massless, sufficient phase space is 
available for the hG and ZG modes for a Higgsino somewhat heavier than h or Z. And this 
mass range will be accessible at the upcoming Run II at the Fermilab Tevatron. The Higgs 
final states also present the exciting possibility of discovering and studying the Higgs boson 
in association with supersymmetry. 

If supersymmetry is broken at a low scale, as required for the di-boson sigatures discussed 
here, it is very likely that the SM gauge interactions play some role in coupling the SUSY 
breaking sector to the SM superpartners 0. However, such gauge-mediated SUSY breaking 
requires additional interactions between the Higgs and SUSY breaking sectors in order to break 
certain Higgs sector global symmetries and obtain acceptable electroweak symmetry breaking 
0. These interactions can modify the Higgsino mass from minimal expectations, and allow 
for a Higgsino NLSP. So searches for di-boson signatures of a Higgsino NLSP within theories 
of low scale gauge-mediated SUSY breaking are very well motivated as possible indirect probes 
for the existence of these additional couplings. 

2 Higgsino decays and production 

The Higgsinos H u and are fermionic superpartners of the Higgs boson fields H u and H^. 
The neutral Higgsinos mix with the gaugino superpartners of the 7 and Z gauge bosons, while 
the charged Higgsino mixes with the gaugino superpartner of the W gauge boson. In the 
limit relevant here, in which the gauginos are heavier than the Higgsinos, the two lightest 
neutralinos and lightest chargino, x?> XsL Xii are predominantly Higgsino and approximately 
degenerate. The splitting between these states is on the order of 10-15 GeV for masses in the 
range 120-250 GeV discussed below. If the U(l)y and SU{2)l gaugino mass parameters, M\ 
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Figure 1: Branching ratios of the lightest neutralino Br(%5 — * G + 7, h, Z) 
as a function of the neutralino mixing angle tan -1 (/x/ Mi), for a fixed mass 
M^o = 160 GeV and m h = 105 GeV for (a) tan/3 = 3 and (b) tan/3 = 40. 



and M2, have the same sign, sgn(MiM2) = + then Xi is the NLSP. For sgn(MiM2) = — it 
is however possible in certain regions of parameter space that xf is the NLSP. In this letter 
only a Xi NLSP, which leads to the interesting di-boson signatures, will be considered. 

The branching ratios Br(x5 ~~ * G + (7, h, Z)) are determined by the Higgsino and gaugino 
content of Xi H> H ■ This is illustrated in Fig. ^ as a function of the neutralino mixing angle 
tan _1 (ii/Mi) for fixed Xi mass, where \x is the Higgsino mass parameter, and tan/3 = v u /vd 
is the ratio of Higgs expectation values. For definiteness the Higgs decoupling limit in which 
decays to the heavy scalar and pseudoscalar Higgs bosons, H and A, are kinematically blocked 
is employed throughout. For gaugino- like Xi the 7 mode dominates, but for Higgsino-like Xi 
the h and Z modes become important. The dependence on sgn(/x) and tan (3 apparent in Fig. [l] 
can be understood in terms of the Xi quantum numbers and couplings and will be presented 
elsewhere. 
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Figure 2: Total cross-section a to t(pp —> XiX^XiXi) i n f° f° r hj = 1)2 times 
the branching ratio into various di-boson final states as a function of the light- 
est Higgsino mass M^o, for a fixed ratio \xjM\ = —3/4, with (a) tan/3 = 3 or 
(b) tan j3 = 40. The center of mass energy is 2 TeV and rrih = 105 GeV. The 
solid line indicates the total cross section. 



The branching ratios also depend on the Xi mass through the phase space available to 
the h and Z modes which suffer a /3 4 velocity suppression near threshold || ||. So even a 
Higgsino-like Xi decays predominantly by Xi ~ * lG for masses not too far above the h and Z 
masses. The mass dependence of the branching ratios is illustrated in Fig. [| in which the pp 
signal cross section times branching ratio into the di-boson final states is given as a function 
of the Xi mass for fixed Higgsino-neutralino mixing. With Xi Higgsino-like the hh, ZZ, or 
hZ modes dominate for very large masses, while the 77 mode dominates for smaller masses. 
However, because of the strong phase space suppression near threshold there is a transition 
region which extends over a significant range of mass between these limits in which the mixed 
final states 'jh and/or 7Z (depending on sgn(/x) and tan/3) are important. These final states 
are particularly useful for masses in the transition region since the photon is quite hard. 

The total cross section cr to t(pp — > XiXji X1X1) f° r hj = 1, 2 in Fig. ^, summed over all the 
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Higgsino-like states, is the relevant signal cross section since these states are approximately 
degenerate, and can all be produced at similar rates. The heavier states cascade decay to Xi 
through neutral and charged current interactions. The partons from these cascade decays are 
relatively soft and probably not particularly useful at the trigger level. 

3 Z boson final states 

The final states with a Z boson can be significant for large tan /5, or at small tan f3 with /x > 0. 
The Z boson can decay invisibly, leptonically, or hadronically, Z —> uu, ££,jj, leading to many 
possible signatures. The ee and fifi leptonic decays allow the possibility of precise reconstruc- 
tion of the Z invariant mass, but suffer from small branching ratio, Br(Z — > ee, /i/i) ~ 6.7%. 
In contrast, the invisible and hadronic decay modes can be useful because of larger branching 
ratios, Br(Z vv) ~ 20%, and Br(Z -» jj) ~ 70%. 

The '-yZIpx di-boson mode dominates the total cross section in the transition region of 
masses as shown in Fig. 2(b). Leptonic decay of the Z provides the cleanest final state, 
, y£ + £~^T, which is similar to existing SM Z'j studies without I^t 10, 0- For a Higgsino search, 
however, an additional large I^t cut, as well as a more stringent photon Et cut should reduce 
the backgrounds to a negligible level. Our Monte Carlo estimates indicate that this channel 
is practically background free, but is limited by the small leptonic branching ratio of the Z 
boson. The Tevatron Run Ha with 2 fb _1 of integrated luminosity will have a reach at the 3a 
discovery level for Xi masses up to 155 GeV for the parameters of Fig. 2(b), while the reach 
in Run lib with 30 fb _1 should approach 220 GeV. 

Invisible decay of the Z gives rise to the signature 7$r- This channel has been studied in 
Run I as a probe for anomalous ■yZ couplings || Backgrounds include jj and jj with one 
jet faking a photon and in each case the remaining jet energy mismeasured to be below the 
minimum pedastool. The largest background in Run I was from single W production with 
W — > ev and the electron misidentified as a photon. This background can be substantially 
reduced by raising the photon Et and J^t cuts above 50 GeV, beyond the Jacobian peak for 
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W — > £v [[H]]. This also reduces the hadronic background. The 3a discovery reach in x\ mass 
should then approach 150 (185) GeV in Run Ha (lib) for the parameters of Fig. |2](b). 

Hadronic decay of the Z in the jZfa mode gives rise to the signature jjjfa. Backgrounds 
are similar to those of the 'jfa channel. The ^yjjfa channel has been studied in Run I in 
order to place limits on squark and gluino masses in very specific supersymmetric models 
|TT| 1. Further background suppressions not included in the Run I study are possible with 
acoplanarity, sphericity and invariant dijet mass cuts to reconstruct the Z boson, and a lepton 
veto. In any case, the total background is expected to be smaller than for the jfa channel, 
due to the presence of two additional hard partons. Given the significant Z hadronic branching 
ratio, the 'jjjfa channel should provide somewhat better reach than the j£ + £~fa or jfa 
channels in Run II. 

The ZZfa di-boson mode dominates at larger x\ mass as shown in Fig. 2(b). Leptonic 
decay of each Z boson gives rise to the spectacular signature £ + £~£' + £'~ fa, with the lepton 
pairs reconstructing the Z mass (in one choice of pairing for £ = £'). This channel is expected 
to be essentially background free, but suffers from small leptonic branching ratio. Because 
of this Run lib will not be sensitive to this channel for the paramters of Fig. |2|. But for 
///Mi = 1/3 and tan/3 = 3 with larger Br(xi —> ZG) (c.f. Fig. [I]), the 3a discovery reach 
in Run lib for the Xi mass is 170 GeV. At the LHC £ + £~£'+£'~fa would represent the gold 
plated channel for the ZZfa di-boson mode from Higgsino decay. 

Hadronic decay of one of the Z bosons gives the signature £ + £~jjfa. An important 
background in this channel comes from ti production with t —>■ Wb and W — ► £v with the 
£ + £~ pair reconstructing the Z mass, and each 6-jet not identified as a heavy flavor. Other 
backgrounds arise from ZZ and W Z in association with jets. In Run lib the 3a discovery reach 
in x\ mass should approach 195 GeV for fi/Mi = 1/3 and tan/3 = 3. Rejecting backgrounds 
for the other decay channels of the ZZfa di-boson mode presents more serious challenges. 
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4 Higgs boson final states 

The decay of Higgsinos to real Higgs bosons gives perhaps the most interesting di-boson final 
states because of the opportunity to study both supersymmetry and the Higgs sector. Higgs 
boson final states are important for small tan /3 and \x < or for large tan /3 with sufficiently 
large Xi mass, as shown in Figs. H and |||. 

In the transition region of Xi mass, jh^x is the most important di-boson mode. With the 
dominant decay h — > bb this leads to the signature jbb^T- Backgrounds include Zjj and Zjj 
with Z — > bb and bbyj and bbjj with one jet misidentified as a photon and in each case the 
remaining jet energy mismeasured to be below the minimum pedastool. Based on the work 



presented here g]] it has been estimated |n| that with a single 6-tag the 3a discovery reach in 
Xi mass should approach 210 (250) GeV in Run Ha (lib) for the parameters of Fig. 0(a). 

For larger x\ masses the hZlpT and/or hh^T modes can become important, as shown in 
Fig. 0. The hhlpx di-boson final state gives rise to the signature bbbblp^- The sizeable QCD 
and electroweak backgrounds to this final state can be significantly reduced by requiring at 
least 3 tagged 6-jets with large invariant mass for two 6-jet pairs [|[], as verified by Monte Carlo 
simulation [13]. Remaining backgrounds include ZZj with each Z — > bb, bbjj with one jet 
misidentified as a 6-jet, and bbbbj with in each case the jet energy mismeasured to be below 
the minimum pedastool, and ti production with t —> Wb and one hadronic decay W — > jj 
with one jet misidentified as a 6-jet, and one leptonic decay W — > £v with £ not identified. 
Accounting for the ti background [13|], the 3cr discovery reach in Xi mass at Run lib should 
approach 240 GeV for the parameters of Fig. ^|(a). 

The hZIpT di-boson mode arising from Higgsino decay is similar to direct hZ production. 
Invisible decay of the Z gives the signature bb$T, and would contribute slightly to searches for 
the SM Higgs boson in this channel. Leptonic decay of the Z gives the signature £ + £~bblfiT- 
Unfortunately, the dominant background from ti production with t Wb and W — > Iv with 
the £ + £~ pair reconstructing the Z mass, is very similar to the signal. Because of this, Run 
II is not expected to be sensitive to this channel. Hadronic decay with Z — > bb gives the 
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Figure 3: Signal cross-section times branching ratio contours in fb for the (a) 
'-fbb^T and (b) bbbb^T channels, as a function of the neutralino mass M^o, and 
the Higgs mass m^, for tan/3 = 3 and fi/Mi = —3/4. 



signature bbbbJ^T, similar to the hhlpx mode. However, because of the smaller branching ratio, 
Br(Z — > 66)/Br(/i — > bb) ~ 20%, Run II will just marginally not be sensitive to the hZI^T 
mode in this channel for the parameters of Fig. |^. 

The Higgs boson final states of Higgsino decay discussed above present the possibility of 
collecting a relatively clean sample of events which contain real Higgs bosons. It is therefore 
interesting to consider the reach as a general function of both Higgsino and Higgs masses. 
The total cross section times branching ratio contours for the ^bblpx and bbbblpT channels as a 
function of the h and Xi masses are shown in Fig. 3. These contours include Br(%5 — > {l, h)G) 
for tan/5 = 3 and ///Ml = —3/4 and SM values for Br(/i — > bb). The Run Ha 3a dicovery 
reach quoted above for the 766^ channel corresponds to a signal times branching ratio cross 
section of 5 fb. For the parameters of Fig. 3 this corresponds to a Higgs mass of up to at 
least 120 GeV for x? masses in the range 135-200 GeV, with a maximum reach in Higgs mass 
of just over 130 GeV. This is to be contrasted with the search for the SM Higgs from direct 
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Wh and Zh production. These SM channels are background limited, and no sensitivity to a 



Higgs mass beyond current limits is expected in Run Ha ||14|| . So the ^bblpT channel presents 
the interesting possibility for Run Ha of a SUSY signal which contains real Higgs bosons. The 
Run lib 3cr dicovery reaches quoted above for the ^bblpT and bbbbl^T channels correspond to 
signal times branching ratio cross sections of 1 fb and 4 fb respectively. For the parameters 
of Fig. 3 the maximum reach in Higgs mass then corresponds to just over 145 GeV and 115 
GeV respectively. 

In order to identify the Higgs boson directly in a sample of events arising from Higgsino 
decays it is necessary to observe a peak in the bb invariant mass. The identifiable di-boson 
final states and large $t carried by the Goldstinos render the supersymmetric Higgs boson 
final states discussed here relatively clean. Reconstructing the Higgs mass peak should be 
relatively straightforward compared to SM Wh and Zh production modes which suffer from 
much larger continuum bb backgrounds. 

All the new signatures presented here involve hard photons, leptons, and/or 6-jets, in 
association with significant missing energy. New triggers are therefore not required, but final 
state specific off line analysis should be implemented in order to search for supersymmetry 
and/or the Higgs boson in these interesting channels. 

Finally, Higgsino decay with a measurable macroscopic decay length to the Goldstino would 
render all the di-boson final states discussed here essentially background free. A search for 
such final states requires a special analysis for displaced £ + £~, jj, or bb with large invariant 



mass and approximately uniform angular distribution with respect to the beam axis fl5 
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